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approximate scaling observed for different hadrons in high energy hadronic 
collisions are discussed. 



Self-Organized Criticality paradigm is a plausible picture for hadron produc- 
tion. A power-law behavior of hadron transverse momentum spectra and an 
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D . Multi-particle production in high energy hadronic collisions is understood 

in Quantum Chromodynamics as a result of irradiation and recombination 
of many strongly interacting quarks and gluons. A system of hadronizing 
^ ■ soft partons is an open dynamical system (the number of partons is not con- 

served) with many degrees of freedom which is far from equilibrium. Even 
though we can calculate exactly an elementary parton-parton interaction, a 
complexity of the system of hadronizing partons doesn't allow the extension 
of these calculations to predict particle states, in which the vast majority 
of the particles are produced with low transverse momenta {Pt)- In addi- 
tion, cascading decays of produced hadronic resonances make the picture 
even more complex. In practice, particle production in high energy hadronic 
collisions is modeled using statistical techniques. 

The seemingly structureless multi-particle final states, however, are found 
to display regularities, namely, long-range particle-particle correlation and 
self-similarity in the particle spectra. Interestingly, a similar behavior is 
found for a broad class of phenomena including turbulence, earthquakes, 
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forest fires, cloud formation, etc. The phenomena from this class are de- 
scribed by power-law distributions. To explain the occurrence of power-laws 
in nature P.Bak et al [|I|] have proposed a Self-Organized Criticality (SOC) 
paradigm. In the SOC the dissipative systems spontaneously evolve into 
barely stable critical states, with long range correlations. The critical state 
is distinguished from the fully deterministic states by the response of a system 
to an external perturbation. A reaction of deterministic system is described 
by a characteristic response at a given scale depending on the system's param- 
eters. The distribution of the resulting events is narrow and is well described 
by an average value. For a critical system, the same perturbation can lead 
to an event of any scale. Although large power events are comparatively 
rare, the mechanism is the same to explain the rare large events and the 
smaller, more common ones. This scaling behavior is described by a power 
law distribution. 

In it's relatively short history SOC patterns have been found in myriads 
of fields including biology, sociology, computer science, economics, etc. The 
SOC models are widely used in modern physics to describe processes, such 
as occurrence of solar flares, earthquakes, and windstorms. The SOC gives 
rise to fractal dimensions of natural objects. The ubiquitous "fingerprints" 
of the SOC states are a power-law distribution of the events and long-range 
correlations. In high energy physics the long-range correlations in particle 
production have been found and studied in details before 0. In the present 
paper we discuss a power-law behavior of the transverse momentum spectra 
of different hadrons measured in collider experiments. 

At low transverse momenta an exponential behavior of the hadronic spec- 
tra has been observed. It was interpreted using a thermodynamic analogy, 
with single particles distributed within a hot hadronizing matter according 
to R. Hagedorn || 



~^m* + P?-exp(- y * ), (1) 

y=0 

where m is a particle mass and T is characteristic temperature of the interac- 
tionQ. The data on high-Pi particle production show a significant deviation 

1 Since the most of the data on hadron production are available for central rapidi- 
ties (yiab ~ 0) we consider the invariant differential cross sections for central rapidity 
only. 
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from the exponential form towards a power-law 



d 3 <7 

E ~ , 

d 3 p \P t 



(2) 



anticipated by the perturbative QCD calculations for parton scattering. These 
calculations, however, fail in the low-P 4 region. An alternative statistical 
hadronization model j|] based on the formation and decay of pre-hadronic 
clusters was found to give a good data description within a broad range of 
values of hadron P t , though it fails to describe the interactions with collision 
energy larger than 30 GeV ||. To describe high energy data within the whole 
P t interval Gazdzicki et al || have suggested a new statistical parameteriza- 
tion of hadron spectra which approximates the exponential form at low-P t 
and the power-law at high-P t 
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In practice, to describe the hadron spectra an empirical power-law func- 
tion is widely used 
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The form (01) is a very close approximation of the following expression 
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(T + (q-l)Jm 2 + P t 2 ) — 



(5) 



obtained in a framework of non-extensive statistical mechanics J7L |8| . In ([5]) 
the T is a temperature of hadronizing matter, m is a hadron mass and q is 
an entropic index. For q — > 1 the expression (|5]) reduces to the Hagedorn 
form ([!]). The non-extensive statistical mechanics generalizes the standard 
thermodynamic approach for a presence of a long-range strong interaction 
and is deeply related to the SOC |§. It was shown JTOl, §] that non-extensive 



statistical model describes the hadron spectra in e + e~ interaction. It is of 
interest to investigate which parameterization from (Q) to (|5]) is preferred by 
the high energy hadron collider data. 

To compare the quality of the data description by the forms (|3|)-(|5|) we 
fitted these forms to the charged particle invariant cross section measured 
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Figure 1: The measured charged particle P t distribution with an overlayed 
the power-law curve (a), and a relative difference between the data and the 



fit results using power-law functions of (Jm 2 + P 2 ) (b) and (P + Pt) (c). 



in the UAl experiment These data are selected as they cover a broad 
interval of Pt values starting from very low P t . The UAl data are shown in 
Fig. 1(a) as function of Pf. Fig lb and Fig lc show the relative differences be- 
tween the data and the fit results using form ([3]) with free mass parameter m 
and the form (^) correspondingly. It is seen from the Fig. 1 that the form (||) 
doesn't describe well the shape of the charged particle P t distribution, while 
the power-law form (f|) is in a good agreement with the data over almost 
whole Pt interval except of the very high P t . Similarly good agreement with 
the data was found for the expression @, however, for simplicity, we use the 
form (H) only in the rest of the paper. 

The power-law behavior of the data is a formal signature of the SOC 
pattern. The SOC model itself doesn't predict the values of the parameters 
Po and D. These parameters are defined by the underlying QCD dynamics. 
The exponent D has to match the pQCD calculations at large Pt- Since QCD 
is flavor blind the exponent D is likely to be the same for different hadrons 
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produced at the same collision conditions. If the SOC approach is applied 
to particle production, P defines a minimal scale at which the P t scaling 
starts to be broken and is likely to be related to the mass (as the only one 
dimensional macroscopic parameter available) of the produced hadron. To 
test this hypothesis we compare the spectra for different long-lived hadrons 
assuming, for simplicity, Pq in to be equal to the hadron mass m. A 
compilation of the invariant cross sections for 7r + [0], K + , A fT3| , |l4j 
and D* + |15j measured in photoproduction at HERA {-Js 1P ~ 200 GeV) is 



shown in Fig.2 as function of (m + Pt). To make the comparison of different 
species of hadrons the cross sections are given for the particle's isospin and 
spin projections. The pions were not identified in the HERA experiments 
and have been recalculated from the measured charged particle spectrum by 
reducing this spectrum by 40% to take into account an admixture of kaons, 
protons and long lived charged leptons. Surprisingly, the spectra shown in 
Fig.2 for different hadrons can be described with a good approximation by 
single power-law function^ The production mechanism for different long 
lived particles and for different intervals of P t is self-similar, i.e. once estab- 
lished for charmed mesons, this mechanism is valid for light pions also at 
low- Pf We test further this universality with the high energy pp interaction. 
In Fig. 3 the invariant cross sections for n + [17]], K° D* + [[0| and B + [20[ 



measured in pp collisions at the Tevatron (y/s = 1800 GeV) are presented 
as function of the (m + P t ) scaling variable. The Tevatron data show an 
approximate scaling behavior for the it, K and D* mesons similar to that 
found at HERA 7p collisions. 

As seen from the Fig. 3 the 5-meson invariant cross section is a factor of 
10 higher than that of the lower mass hadrons when plotted at the same val- 
ues of (m + P t ) variable. At HERA, though the B mesons are not yet directly 
reconstructed, the 5-meson inclusive cross section was estimated using the 
Lund model [^TJ with a normalization to the 6-quark production cross section 



measured at HERA [ 22"| . The result of these calculations is shown in Fig.2 
with a shaded band. The band width is defined by the experimental uncer- 
tainty of the 6-quark production cross section measurement. The calculated 
fi-meson yield at HERA is again a factor of 10 higher than that expected for 
other hadrons shown in Fig. 2. Noteworthy, the QCD calculations for open 
beauty production strongly underestimate the measured cross section at the 



2 Recently, the HI has reported ||16| the preliminary proton photoproduction cross sec- 
tion which supports the observed scaling. 
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Tevatron and HERA as well. Recently, a new parton ^-factorization 
approach was used to describe the open beauty production p4| , but it is 



still to be demonstrated that this approach can be concistently applied to 
other hadron spieces. Within the statistical approach the observed excess of 
.B-meson production over the universal behaviour of the scale dependence of 
other mesons signals a different dynamics or the existence of an additional 
mechanism of 5-production. 

In summary, we noted that particle production in high energy hadronic 
collisions has the formal properties of a Self-Organized Critical process, name- 
ly: a) the hadronizing partons are known to represent an open strongly inter- 
acting dynamical system far from equilibrium - a typical system to deal with 
the SOC; b) the resulting hadronic state shows long-range particle-particle 
correlations; c) the particle spectrum is described by a power-law distribu- 
tion as function of transverse momentum and mass of the particle; d) this 
distribution is the same for different long-lived hadrons and is defined by 
macroscopic parameters like mass and spin of the particle as expected for 
statistical models. Whether the latter is a pure coincidence or a demon- 
stration of a statistical (alike SOC) nature of the hadron production is an 
open question. To prove an eligibility of the SOC mechanism in particle pro- 
duction, a successful quantum cellular automaton describing a quark-gluon 
system has to be built. It is possibly could be done in the future using QCD 
lattice calculations. 

For completeness we give also the earlier applications of the SOC in high 
energy physics known to us. The SOC was suggested |25| to describe inelastic 
diffractive scattering. Fractal properties of the hadronic collisions have been 
discussed in ]26[ . In |27[] a concept of fractal dimensions of the proton struc- 



ture function was introduced. This gave an excellent description of the low-x 
HERA data, both in the non-perturbative and the deep-inelastic domain. 
These observations together with the self-similarity in particle production 
discussed in this paper make the SOC approach to be a candidate for an 
effective picture to describe a variety of non-perturbative QCD phenomena. 
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Figure 2: The invariant cross- sections for different long-lived hadrons as 
function of m + px, for rapidity y = as measured in the laboratory system 
at HI and ZEUS detectors. The cross-sections are given for one spin and 
isospin projections. The shaded bar corresponds to the S-meson cross section 
estimated from the measurement of inclusive 66-production. 
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Figure 3: The invariant cross-sections for different long-lived hadrons as a 
function of m + pt, for rapidity y = as measured in the laboratory system 
at CDF detector. 
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